The Anderson transition in quantum chaos 
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We investigate the effect of classical singularities in the quantum properties of non-random Hamil- 
tonians. We present explicit results for the case of a kicked rotator with a non-analytical potential 
though extensions to higher dimensionality or conservative systems are straightforward. It is shown 
that classical singularities produce anomalous diffusion in the classical phase space. Quantum me- 
chanically, the eigenstates of the evolution operator are power-law localized with an exponent uni- 
versally given by the type of classical singularity. For logarithmic singularities, the classical motion 
presents 1// noise and the quantum properties resemble those of an Anderson transition, namely, 
multifractal eigenstates and critical statistics. Neither the classical nor the quantum properties de- 
pend on the details of the potential but only on the type of singularity. It is thus possible to define 
new universality class in quantum chaos by the relation between classical singularities (anomalous 
diffusion) and quantum power-law localization. 

PACS numbers: 72.15.Rn, 71.30.+h, 05.45. Df, 05.40.-a 
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The description of quantum features based on the un- 
derlying classical dynamics has been intensively investi- 
gated since the early days of quantum mechanics. For 
instance, the celebrated Bohigas-Giannoni-Schmit con- 
jecture U states that the spectral correlations of a quan- 
tum system whose classical counterpart is fully chaotic 
depend only on the global symmetries of the system and 
are identical to those of a random matrix with the same 
symmetry (usually referred to as Wigner-Dyson statis- 
tics (WD)). However this conjecture is not always ver- 
ified. A paradigmatic example is the kicked rotor: Al- 
though classically chaotic, its spectral correlations follow 
Poisson statistics and eigenfunctions are exponentially 
localized. This behavior was fully understood after 
mapping the kicked rotator problem onto an ID disor- 
dered system (free particle in a random potential). It 
is by now well established that, in less than three di- 
mensions, all the eigenstates of a disordered system are 
exponentially localized in the thermodynamic limit for 
any amount of disorder. However, in three and higher di- 
mensions, there exists a metal insulator transition (MIT), 
also called Anderson transition, for a critical amount of 
disorder. At the MIT, the wavefunction moments V q 
present anomalous scaling with respect to the sample 
size L, V q = J d d r\vj(r)\ 2 i cx L-- ^?- 1 ) , where D q is 
a set of exponents describing the transition. Wavcfunc- 
tions with such a nontrivial scaling are said to be mul- 
tifractal 0]. Spectral fluctuations at the MIT are com- 
monly referred to as 'critical statistics' yL Typical fea- 
tures include: scale invariant spectrum [5j, level repul- 
sion and linear number variance [|J. Similar properties 
have also been observed in disordered systems with ran- 
dom power-law hopping (linked to classical anomalous 
diffusion ||) provided that the exponent of the hopping 
decay matches the dimension of the space. In the ID case 



(1/r decay) the relation to a MIT has been analyti- 
cally established. For non-random power-law hopping a 
MIT was recently reported 01 even for power-law expo- 
nents larger than the dimensionality of the space. Unlike 
the random case, critical states appear only in a certain 
energy window. A similar power-law behavior has also 
been found in certain non-random Hamiltonians: Fermi 
accelerator fill. Coulomb billiard 0, Anisotro pic Ke- 
pler problem |l3| and generalized Kicked rotors [lil Hif . 
In all of them the potential has a singularity and con- 
sequently the KAM theorem does not hold. The clas- 
sical dynamics of these systems is intermediate between 
chaotic and integrable and the classical phase space is 
homogeneously filled with cantori. The goal of this letter 
is to show that, in certain cases, classical chaotic Hamil- 
tonians with non-analytical potentials present classical 
anomalous diffusion and quantum power-law localization 
of the eigenstates. The case of logarithmic singularities, 
corresponding to a MIT, will be investigated in detail. 
Classical mechanics-We study the Hamiltonian 



H = P - + V{q)Y,5{t-nT) 



(1) 



with V(q) = e\q\ a and V(q) = elog(M); q € [-tt.tt), 
a G [— 1, 1] and e a real number. The dynamical evolu- 
tion over T is dictated by the map: p n +i = Pn ~ 9 V g q qn ^ i 
Qn+i = q-n + Tp n+ i (mod 2tt). Our aim is to show that, 
though obviously deterministic, the motion is accurately 
modeled by an anomalous diffusion process. In particular 
we compute the classical density of probability P(p, t) for 
finding the particle at time t with a momentum p. Un- 
like the case of normal diffusion, the information obtained 
from the knowledge of a few moments of the distribution 
is not sufficient to fully characterize the classical motion. 
For instance, the second moment may be (r 2 ) ~ t but 
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FIG. 1: P(p,t) versus p (log scale). In all cases t = 2000 
except for a — —0.5 where t = 40. For V(q) = \q\ a with 
a = 0.9 (solid) and V(q) = (\q\ + 0.005) 04 (dash) the diffusion 
is normal (squares and circles are for best Gaussian fitting 
respectively). However for a = 0.4 (dot) and a — —0.5 (dash- 
dot), P(p,t) oc tp~'( a \ Inset shows 7 as a function of a. The 
star corresponds to V(q) = log ( | g | ) (see Fig. 2). 
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FIG. 2: P(p,t = 40) versus p (semi-log scale). For V(q) = 
log(|g|) (solid) P(p,t) is very well fitted by the Lorentzian 
distribution P(p,t) — 7 i+pffs m with v(t) = 2/t (squares). 
As shown in the inset, P(p,t) (solid) is not affected by adding 
a noise term cos(q) to V(q) (circles). 



this by no means assures that the density of probability 
is Gaussian like (lif . 

The classical density of probability P(p, t) has been 
evaluated by evolving 10 s initial conditions uniformly dis- 
tributed in (q,p) G [— ir, n) x [— tt, 7t). For a > 0.5 (see 
Fig. 1), the diffusion is normal and P(p,t) ~ e~ cp '* l\fi 
(c a constant). However, for —0.5 < a < 0.5 and t > p, 
P(p, t) ~ tp 1 ^ has power-law tails (see inset in Fig. 
1 for a relation between 7 and a). Such tails are con- 
sidered a signature of anomalous diffusion. They are 
indeed a solution of the fractional diffusion equation, 
(I " iS^r) P iP^) = S(p)S(t) (see Q for the def- 
inition of fractional derivative). In order to explicitly 
show that the origin of the anomalous diffusion is ex- 



clusively related to the classical singularity we have re- 
peated the calculation with the singularity-free potential 
V(q) = e(\q\ + b) a , b > 0. As expected, the diffusion 
goes back to normal as in the kicked rotor with a smooth 
potential (see Fig. 1). We have also checked that our 
results are generic, namely, they do not depend on the 
details of the potential but only on the kind of singu- 
larity. Thus P(p, t) is not altered by adding a smooth 
perturbation V per in Eq.QJ (see inset of Fig. 2). The po- 
tential V(q) — elog(|q|) has been studied in detail since 
its quantum counterpart (see below) resembles a disor- 
dered conductor at the MIT. In this case (see Fig. 2) 
P(p, t) is accurately fitted by a Lorentzian. We recall this 
special form of P(v , t) describes the dynamics of systems 
with 1/f noise 17]. In passing we mention that in this 
case the classical phase space is dominated exclusively 
by cantori with no coexistence of chaotic and integrable 
regions. 

Quantum mechanics- We now discuss the quantum prop- 
erties of Eq. Q . For the sake of clearness let us first state 
our main results: 

1. The eigenstates of the evolution matrix associated to 
the Hamiltonian Eq. are power-law localized with a 
exponent depending only on the singularity a. 

2. For V(q) — elog(|g|), the eigenstates are multifrac- 
tals and the spectral correlations are described by critical 
statistics as at the MIT. 

The quantum dynamics of a periodically kicked system 
is governed by the quantum evolution operator U over 
a period T. Thus, after a period T, an initial state tpo 

-ip 2 T iV(i) -ip 2 T 

evolves to ^(T) = Ui/jq = e e ^ e *>> ipo where 
p and q stand for the usual momentum and position 
operator. Our aim is to solve the eigenvalue problem 
U^n = er ltin / h ^i n where ^ n is an eigenstate of U with 
eigenvalue K n . In order to proceed we first express the 
evolution operator in a matrix form (m\U\n) in the basis 
of the momentum eigenstates = ^75^ }• For practical 
calculations one option is to make the Hilbert space finite 
(i.e. m, n = 1, . . . N) but still keep the matrix Unitary. 
This can be achieved by imposing a periodic condition 
(with period P) to the momentum. A large period makes 
the phase space effectively resemble a cylinder (P = N 
in this paper; for more details see 0])- I n or der to keep 
the kinetic term of the evolution matrix also periodic, we 
take T = 2irM/N with M an integer (not a divisor of N). 
M is specified to make T roughly constant (w 0.1 in this 
paper) for every N used. Based on the relation between 
the period in momentum P and the period in position Q 
(PQ = 2ttNTi), h = 1. The resulting evolution matrix 
(for N odd) then reads 

(m\U\n) = L e -^Mn 2 /N ^ e m,m,n) (2) 

where 4>{l,m,n) = 2ir(l + 9 )(m - n)/N - iV(2ir(l + 
9 )/N), l = -(N - l)/2, . . . (N - l)/2 and < O < 1; 
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for systems with broken time reversal invariance, 
n 2 h 2 sin 2 (7rs) 



R 2 {s) = -K\a) 



FIG. 3: The level spacing distribution P(s) in units of the 
mean level spacing s = (Ei — i?i_i)/A. A transition is ob- 
served from Poisson to WD statistics (keeping e constant) as 
a goes from negative to positive. In both cases N = 3100. 




FIG. 4: The spectral rigidity As(L) versus L. Symbols are the 
numerical calculation for V(q) = elog(g). Lines are the ana- 
lytical prediction Eq.JSJ. The right (left) panel corresponds 
to the case of (broken) time reversal invariance. 



0q is a parameter depending on the boundary conditions 
(0q = for periodic boundary conditions). The eigenval- 
ues and eigenvectors of hi can now be computed by using 
standard diagonalization techniques. For 9q = 0, parity 
is a good quantum number and consequently states with 
different parity must be treated separately. 

We first study the spectral fluctuations of the unfolded 
spectrum. The analysis of short range spectral correla- 
tions as the level spacing distribution P(s) indicates (see 
Fig. 3) that for a > 1/2 (a < -1/2), P(s) tends to Pois- 
son (WD) as the volume is increased. In the intermedi- 
ate region —1/2 < a < 1/2 long range correlators as the 

spectral rigidity A 3 (L) = Jq(L 3 — 2L 2 x + x 3 )Y< 2 (x)dx 

where £ 2 (L) = (L 2 ) — (L) 2 = L + 2 (L — s)R 2 (s)ds 
is the number variance and R 2 (s) is the two level spec- 
tral correlations function, are more suitable for a qual- 
itative analysis due to the possibility of performing a 
spectral average. We restrict ourselves to the potential 
V(q) = elog(g) related to classical 1/f noise. The R 2 (s) 
associated to this type of motion has been evaluated in 



sinh 2 (7r 2 /is/2) 



(3) 



where h -C 1 is related to e by h = 1/e -C 1. In Fig. 
4 we show A 3 (L) for different e's and N = 3099. The 
time reversal invariance in Eq.|J2J is broken by setting 
T = 2tt(3 with [3 an irrational number. As observed in 
Fig. 4 (left), A 3 (L) is asymptotically linear, this feature 
is typical of a disordered conductor at the MIT. More- 
over, the agreement with the analytical prediction based 
on Eq. © is excellent. We remark that the value of h 
best fitting the numerical result is within 5% of the ana- 
lytical estimate h = 1/e. We have repeated the calcula- 
tion keeping the time reversal invariance (T = 2nM/N) 
of Eq.© (Fig. 4, right). We do not have in this case 
an analytical result to compare with, but A 3 (L) is also 
asymptotically linear as in the previous case. The eigen- 
states >F n = i>n{k)\k) of the evolution operator U are 
also strongly affected by the classical singularity a. In- 
stead of the exponential localization typical of the stan- 
dard kicked rotor, our numerical results indicate that 
1^(^)1 ~ \k — fc m aa:| _( - Q+1 ' ) decays with power- law tails 
(kmax is where |?/> ra (fc)| takes the maximum). This is in 
agreement with the findings of [||. The eigenvector cal- 
culation was carried out in the cylindrical phase space 
(i.e. we switch off the periodic condition in momentum). 
Fig. 5 shows the results after a proper average over all 
N eigenstates. Remarkably, both classical (anomalous) 
dynamics and quantum (power-law) localization are con- 
trolled by the classical singularity a. This fact allows 
us to define a new universality class in quantum chaos 
labeled by a. For V(q) = elog(|g|), \tp n (k)\ ~ 1/fc and, 
according to [§], the system is at the MIT. We now inves- 
tigate how transport properties are affected by the non- 
analytical potential. We compute the quantum density of 
probability P q {k, t) = | (k\<fi(t)) \ 2 of finding a particle with 
momentum kfi at time t for a given initial state |</>(0)). 
We restrict ourselves to the case V{q) — elog(|g|) in or- 
der to compare our findings with the ones at the Ander- 
son transition where it is reported ,19j that, in the limit 
k «: t > 1, P q (k,t) ~ t~ D ' 2 kP' 2 ~ 1 where D 2 is the mul- 
tifractal dimension defined previously. As expected, we 
have observed a similar power-law behavior in our model 
(see Fig. 6). The numerical value of D 2 increases as e 
does and it is within 10% of the one obtained by the in- 
verse participation ratio (not shown) . These results show 
that quantum mechanically the diffusion is also anoma- 
lous but, unlike the classical case, the decay of P q (k,t) 
depends not only on the classical singularity a but also 
on the coupling constant e. The overall effect of the quan- 
tum corrections is also to suppress the classical (anoma- 
lous) diffusion at a rate that increases as e gets smaller. 
To summarize, quantum 1+1D systems with classical log 
singularities share all the properties of a disordered con- 
ductor at the MIT: Multifractality, critical statistics and 
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FIG. 5: Averaged modulus of the eigenstates of the evolu- 
tion operator for N = 2100. We set e = 0.01 to a small 
value in order to avoid the finite size effects. The best fit- 
ting power exponents —0.49, —0.98 and —1.36 correspond to 
a — —0.5, 0(log), 0.4 respectively are very close to the analyt- 
ical prediction — (1 + a) [9j. 
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FIG. 6: Time averaged P g (k,t) over 4950 < t < 5050 and 
P q (k = 0,t) versus t (inset) for V(q) = 101og(|<j|). The best 
fitting (dotted line) is of a slope —0.34 and —0.67 (inset) re- 
spectively, which correspond to a D2 ~ 0.67. 



quantum anomalous diffusion. Analytical results can in 
principle be obtained by mapping Eq.Q onto an ID An- 
derson model. This method was introduced in [2j for the 
case of a kicked rotor with a smooth potential. We 
do not repeat here the details of the calculation but just 
state how the ID Anderson model is modified by the 
non-analytical potential. It turns out that classical non- 
analyticity induces long-range disorder in the associated 
ID Anderson model, T m u m + J2 r U r u m +r = Eu m where 
U r ~ A r /r 1+a is the Fourier coefficient of — ta,nV(q)/2, 
T m = 6 /2 - Tm 2 /2 and E = -U . For random Gaus- 
sian A r this model corresponds to the one studied in 
which is solved by using the supersymmetry method. For 



the special case of a constant non-random A r (closer to 
our case) it was found that 0] some of the eigenvectors 
remain critical even for 1 < a. This suggests that the 
results reported above for a = (log potential) may be 
extended to other cv's though clearly more work is needed 
to further clarify this issue. In conclusion, we have shown 
that classical singularities in chaotic Hamiltonians may 
lead to classical anomalous diffusion and quantum power- 
law localization. Both quantum and classical features are 
governed by the classical singularity. We thus suggest a 
new universality class in quantum chaos labeled by the 
type of classical singularity. In the case of logarithmic 
singularities the classical dynamics presents 1// noise. 
Quantum mechanically the system posses all the features 
of a MIT: multifractal wavefunctions, critical statistics 
and quantum anomalous diffusion. 
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